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Abstract 

We present a general theory of the corrections to the asymptotic behaviour of the 
Renyi entropies Sj? = (1 — n) _1 logTrp^ which measure the entanglement of an 
interval A of length I with the rest of an infinite one-dimensional system, in the 
case when this is described by a conformal field theory of central charge c. These 
can be due to bulk irrelevant operators of scaling dimension x > 2, in which case 
the leading corrections are of the expected form £~ 2 ( x ^ 2 ) for values of n close to 1. 
However for n > x/(x — 2) corrections of the form l 2 - x ~ x / n anc l g-2x/n ar i se anc [ 
dominate the conventional terms. We also point out that the last type of corrections 
can also occur with x less than 2. They arise from relevant operators induced by the 
conical space-time singularities necessary to describe the reduced density matrix. 
These agree with recent analytic and numerical results for quantum spin chains. 
We also compute the effect of marginally irrelevant bulk operators, which give a 
correction 0((log^)~ 2 ), with a universal amplitude. We present analogous results 
for the case when the interval lies at the end of a semi-infinite system. 



1. Introduction 



Recently there has been much interest in characterising bipartite quantum entanglement 
of pure states in extended systems near a quantum critical point in terms of the Renyi 
entropies[T]. For a given division of the Hilbert space into a part A and its complement 
B, these are defined as 

= (l-n)- 1 \ogTrp n A , 

where pa = Tr# p is the reduced density matrix of the subsystem A, and p = \^/)(^/\ is 
the density matrix of the whole system in a pure state \^f). Knowledge of the for 
different n characterizes the full spectrum of non-zero eigenvalues of pA (see e.g. [2j), and 
gives more information about the entanglement than the widely studied von Neumann 
entropy S^. It also gives a fundamental insight into understanding the convergence and 
scaling of algorithms based on matrix product states [3]. 

In [H El E] it was shown that for a one-dimensional critical system whose scaling limit 
is described by a conformal field theory (CFT) of conformal anomaly number (central 
charge) c, in the case where A is an interval of length i embedded in an infinite system, 
the asymptotic behaviour of the Renyi entropies is given by 

fi^la+n-^iogi+oa). (i) 

o 

This result has by now been verified analytically and numerically (see e.g. [7J [8j [9j [TOj 
ITT| [P2| [T3| [TU [15j [16j [T7] , but this list is far from being exhaustive) for a large number 
of examples of quantum spin chains whose scaling limit is believed to be described by 
CFT. It gives, in fact, one of the most accurate ways of measuring the conformal anomaly 
number. However, this asymptotic result is often obscured by large (sometimes oscillating) 
corrections to scaling [TU [191 120] whose origin has, so far, not been clarified in the context 
of quantum field theory. 

In any real system with an ultraviolet cut-off such as a lattice, even when it is tuned 
to the critical point there will in general be operators in the hamiltonian or action which 
ensure that the continuum field theory results are only asymptotic on distance scales much 
larger than the cut-off. It is important to understand the form of the corrections to the 
asymptotic results in order to make accurate fits to analytic and numerical finite-lattice 
data. Renormalization group (RG) theory shows that these are ordinarily due to irrelevant 
operators, with scaling dimension x > d (d being the dimension of space-time), which are 
allowed in the effective action even at the critical point. These are generally responsible for 
finite-size corrections in integer powers 

o{ L -( x -d) j2jj j where L 

is a characteristic length 

scale of the system. Since field theory classifies these irrelevant operators and allows the 
computation of their scaling dimensions, it therefore is capable of quantifying the form 
of corrections to scaling, if not the values of the non-universal prefactors. Moreover, as 
x — > 2+ these corrections become more important, and, at x = 2, the so-called marginal 
case, take the form of logarithms rather than powers. The theoretical advantage is that 
the amplitudes are then often universal and calculable. 

In this paper we analyse the corrections to the Renyi entropies due to irrelevant and 
marginally irrelevant bulk operators in the hamiltonian. For the case of an interval of 
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length i in an infinite system, we find the expected terms 0{t~ 2 ^ x ~^\ For the case of an 
interval at the end of a semi-infinite system there can also be terms 0{t~( x ~ 2 >). However 
we also find unusual terms which are 0{l 2 - x - x ' n ) and 0(£- 2x/n ), in which the exponent 
depends on n. For n > n c (x) = x/(x — 2) these dominate the conventional corrections. 

While the appearance of n-dependent exponents for an irrelevant bulk operator is 
perhaps surprising, in recent analytic and numerical studies of several spin chains [TBj even 
more dramatic corrections of the form £~ 2x / n ^ with values of x < 2, have been reported. 
We argue that these may also be understood from the field theory under the assumption 
that the conical singularities of the Riemann surfaces on which the path integral must be 
evaluated in fact locally break the criticality of the system, thus allowing the appearance 
of operators at these points, which, in the bulk, would be relevant with scaling dimension 
x < 2. Such operators do not drive the system away from bulk criticality because they 
are localised near points in space-time. It is an old result of CFT |22^l that such a bulk 
operator near a conical singularity of degree n in fact has its scaling dimension modified 
to x/n. Thus we predict that such operators should give rise to corrections of the form 
£-2x/n j n cage Q f an j n t erva i j n an infinite system and i~ x l n in the semi-infinite case. 

The marginal case x = 2 turns out to be even more subtle and difficult. However, we 
are able to show that the corrections take a universal form in which c in (JT]) is replaced 
by 

where b is a universal (and, for many CFTs, known) operator product expansion (OPE) 
coefficient. Note that, for this leading correction, the n-dependence is the same as that 
of the leading term. As we discuss below, this is a consequence of Zamolodchikov's c- 
theorem[23j. 

2. Field theory for corrections to scaling 

As was argued in Refs. [U |5j [21] , for a one-dimensional quantum system on a lattice £ the 
expression Tr p n A is given by the ratio Z n /Z™, where Z n is the partition function defined 
by the path integral on n copies of C ® imaginary time r G {— oo, +oo}, sewn together 
along r = so that in the interval B C C the jth sheet with r > is connected to the 
jth sheet with r < 0, while in the interval A the j'th sheet is connected to the (j + l)th 
sheet, cyclically. In the continuum limit, where the lattice is replaced by the real line, 
this defines an n-sheeted Riemann surface TZ n with conical singularities at the ends of the 
interval A, but in what follows later it is important to realise that the identification also 
holds on a spatial lattice. It is also convenient to consider the logarithm —(F n — nF{) of 
the above ratio of partition functions, which is then proportional to the Renyi entropy 

q(n) 
A ■ 

In analysing such a lattice system using continuum field theory methods, the leading 
universal effects of the lattice are usually assumed to be taken into account by imposing 

1 In this reference, only corners on the boundary were in fact considered. However, the analysis extends 
trivially to the case of conical singularities in the bulk. 
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a short- distance cut-off e on the continuum theory: for example by restricting the asymp- 
totic expressions for correlation functions to be valid only for separations greater than e. 
While this might appear to be a rather crude approximation to the effect of a lattice, since 
we are looking only for the universal form of corrections to scaling and not their precise 
amplitudes, this is in fact adequate. The response of the free energy F to a change in e 
is in general given in the field theory by the integrated trace of the stress tensor [25] 126] : 



In a CFT in flat space-time, (Q{z)) = 0, implying that F is scale- invariant. (This is in 
fact strictly correct only after subtracting off the non-universal bulk free energy, which in 
our case automatically cancels in the combination F n — nF%.) However, in Refs. [251 12S] it 
was shown that this is no longer true at conical singularities: in fact each one contributes 
a term (c/12)(n— n ) to the right hand side of ([3]). Integrating up and using the fact that 
F n — nF\ can depend only on the ratio £/e then gives the result (0Q), which was derived in 
slightly different ways in Refs. (3] and [5]. 

Let us now consider the effect of a bulk irrelevant operator in the hamiltonian of 
action. In the field theory this is equivalent to perturbing the CFT by an operator Q(z) 
of scaling dimension x > 2, so the action is 

S = S* + A J $(z)d 2 z, (4) 

where S* is the CFT action, and A is a coupling constant. It has the dimensional form 
g/e 2 ~ x , where g is dimensionless. The change in the dimensionless free energy is then 
formally given by the perturbative series 

- 6F - = E nvT" / " / {Hzi) ■ ■ ■ ^ ZN ^ d " zi ■ ■ ■ d " ZN ' ( 5 ) 

of integrals of connected correlation functions of the CFT over TZ n . 

In the case when A is the interval (0, £) in an infinite system, TZ n may be conformally 
mapped to the punctured complex plane C = C \ {0} by[5j 

/ \l/n C n 

C=( ,) = z = £f(C)=£-^ 



z-£J J ^ J C n -1 

This maps the ends of the interval to £ = and oo. The correlation functions in the two 
geometries are related by 

TV 

. . . Hz N )) nn = J] \tf %){-*(*(&) • • • $((n))c> ■ 
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Since ($)c — 0, the N = 1 term in ([5]) is absent in this case. The second order term may 
be transformed to an integral over the £-plane: 



w n 




l/'(Ci)M/'(C2)l 
IG-C2I 2 * 



2-x 



d 2 (id 2 ( 2 



(6) 



| CiC2 | M( „-i) 



d 2 Cid 2 



(2 • (7) 



c ic icr - 1 



2x ^n _ ^-2x 



ICi - C2 



2.c 



Note that this integral makes sense also for non-integer n, although we derived it only for 
integer n. We may then consider n as an arbitrary real parameter. The integral has now 
several potential sources of UV divergence which should be regulated by a cutoff that in 
the z- plane is 0(e). These potential divergences give rise to a further powers of e in 5Fn 
and, by scaling, to further powers in I in the corrections to scaling. 

To elucidate the mechanism let us start from the case when n — 1 is small and positive. 
The integral in (J7j) then converges for x < 1. For larger values of x, a divergence comes 
from the region d — > ( 2 and it should be regularised with a cut-off \z%— z 2 \ < e. The leading 
divergence in the integral is 0(e 2 ~ 2x ), leading to a dependence in 5Fn^ proportional to 
e _2 Area(7?.„). This is a contribution to the non-universal bulk free energy, which cancels 
in the combination F n — nFi. This subtraction also cancels the apparent singularities at 
( n = 1, which correspond to \z\ — > oo on lZ n . If the leading divergence is subtracted off 
from (J7|), the remainder is analytic at x — 1. In fact, the finite part is then given by the 
analytic continuation of ([7]) around its pole at x — 1. This is then finite all the way up 
to x = 3, and the £ dependence comes from the explicit prefactor, that is £ _2 ( a; - 2 ) ) the 
standard power law of FSS in the RG framework. Note that although the amplitudes of 
these terms are non- universal, depending on g, their ratios are universal, and in principle 
calculable by evaluating the analytic continuation of the integral. We also note that for 
an interval near the end of a semi-infinite system, depending on the boundary conditions 
($) may be non- vanishing, in which case the leading correction will be 0(£ 2 ~ x ). 

However, for larger values of n and x > 2, (J7J) may also exhibit divergences as (j — > 
or oo. In the original coordinates, these occur as z\ or z 2 approach one of the branch 
points. These are genuine divergences due to the conical singularities and are not present 
in the bulk. Close to (say) z = each integral behaves like J \z\^ 1 ^ n ' > ~ 1 ^ x d 2 z, which 
diverges if n > n c (x) — x/(x — 2). In this case the integral must be further regulated 
with a cutoff \z\ < e. This leads to an further multiplicative factor oc e 2 -' x +( x / n ) which, 
once again, by scaling leads to corrections to scaling in the Renyi entropy of the form 
£i-2x-(2-x+x/n) _ £2-x-x/n^ p Qr n > ^ these divergences are further enhanced when z\ 
and z 2 are close to different branch points, leading to a further factor of e 2 ~ x + x / n ^ an d a 
consequent ^-dependence of the form £^-2x-2(2-x+x/n) _ £-2x/n^ ^j^g singularities when z\ 
and z 2 approach the same branch point are removed by the bulk free energy subtraction.) 

To summarise, the presence of bulk irrelevant operators with x > 2 leads to corrections 
to scaling of the form £ 4_2:E , £ 2 - x - x / n an d £~ 2x / n . In general, all these will be present, but 
for n < n c = x/(x — 2) the first dominates, while for n > n c it is the last. For n w n c we 
expect them all to play a role, with multiplicative logarithmic factors when n = n c . 
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The appearance of these terms is similar to what happens in the case of a boundary. 
In that case the bulk operator $(y) at a distance y from the boundary behaves like [27] 
y~ x+Xb &b, where &b is a boundary operator with scaling dimension Xb- This gives rise, on 
integration, to a term \ 2 e X+Xb j n the free energy per unit length of the boundary which 
is singular if Xb < x — 1, and hence to a term L x ~ Xb L 4 ~ 2x in the total free energy of a finite 
system of size L. This boundary contribution can actually overwhelm the normal bulk 
term. Analogously, close to a conical singularity, &(z) ~ \z\~ x+x ' n & n \0), where $ ( - n ) is 
localised at the tip of the cone, and has scaling dimension x/n. The difference is that in 
the boundary case Xb is usually larger than the bulk scaling scaling dimension x, while at 
a conical singularity the scaling dimension (x/n) becomes arbitrarily small for large n. 

3. Relevant operators at conical singularities 

We now argue that corrections to scaling of the form £~ 2x / n in the Renyi entropy can arise 
not only by the mechanism discussed above due to irrelevant bulk operators with x > 2, 
but in a different way which also gives corrections of this form but with x < 2. This is 
because a lattice model which is critical can nevertheless generate operators localised at 
the conical singularities with scaling dimension (x/n) but with x < 2, that is, operators 
which, if they appeared in the bulk hamiltonian, would be relevant and therefore drive the 
system away from criticality. This most easily seen if we discretise time as well as space. 
Consider, for example, a 2d classical model on a square lattice. The anisotropic limit of 
this in general gives rise to a Id quantum hamiltonian. Computing Tr p n A for this quantum 
model corresponds, as before, to considering the 2d lattice model on an n-sheeted surface 
with branch points of degree n. On the lattice, the details of this depend on exactly how 
the degrees of freedom are divided between A and B. In most models like quantum spin 
chains, the degrees of freedom are on the lattice sites, so the division is along a spatial 
bond. In the time-discretised picture, the branch point can be considered to lie halfway 
along a spatial bond, or in the middle of a plaquette. In either case, it is clear that degrees 
of freedom close to the branch points have an enhanced (if n > 1) number of neighbours. If 
we maintain the same bond interactions as in the bulk, there will therefore be an effective 
local coupling to the energy density, or to other operators which do not break the internal 
symmetries of the system. For example, in the case of the nearest-neighbour Ising model, 
if we put the branch point in the centre of a plaquette, each order variable (Ising spin) 
still has exactly 4 nearest neighbours, but the dual spin, which sits on top of the branch 
point, will have 4n nearest neighbours, thus locally breaking the self-duality of the model 
and, locally, driving it away from criticality. In this case we would therefore expect a 
coupling to the energy density, with x = 1. In the field theory, such an operator close to 
a conical singularity of degree n has its scaling dimension modified to x/n |22j, leading, 
in this case, to corrections of the form i~ 2 / n . 

The conclusion is that the correct form of the field theory action on the n-sheeted 
surface should be 

S = S CFT + I <$> 3 d 2 z + x ^l n \p) , 

j p k 
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where the second term takes account of bulk irrelevant operators $j with Xj > 2, already 
discussed, and the third term is a sum over the branch points P of localised operators 
$j^(P) with scaling dimension Xk/n, with all possible values of Xk allowed by symmetry, 
including those with Xk < 2. However, in the perturbative expansion in powers of the 
A fe , each at a given branch point should appear at most once (otherwise we can 

use the OPE to write higher powers in terms of other localised operators). In the case of 
an infinite system, since ($( n )(P)) = 0, we therefore expect the leading correction to be 
0(l~ 2x ^ n ). This should be the case no matter how many branch points > 2 there are. For 
an interval at the end of a semi-infinite system, depending on the boundary conditions it 
may be that (& n \P)) ^ 0, in which case the leading correction will be Oil Xk / n ). 



4. The marginal case and the c-theorem 

Next we turn to the marginal case x — > 2+, which is technically more challenging. The 
integral in ([6]) is in general very difficult to manipulate into a form in which the necessary 
analytic continuation around the pole at x — 1 can be performed. It is easier to extract 
the limit x — > 2 by making the subtraction explicitly, and this we now describe. It is, 
however, quite subtle, as the cut-off and subtraction must be performed in the z-plane. 
Imposing a cut-off \z\ — z%\ > e, in the £-plane we have 

SF n = -h\llef~^ [ [ l/XCi)| 2 -1f(C 2 )| 2 -^ 2Cirf2c2 _ (8) 



/(C0-/(Ca)|>(e/fl 



ICi - c 2 



The subtraction is 



9 2-2x f 9 2-2x 



Tin 



'l<a-Ci|>(e//|/'(Ci)D 




It can be shown that, up to terms which vanish as x — > 2, the cut-off in the last integral 
can be replaced by |/(Ci) — fiCz)] > ( e /^)- Thus the second order contribution to the 
subtracted free energy is 



1/(£/e) 4-2, / / lf(Cl)| 2 -|f(C2)| 2 --|f(Cl)| 4 - 2 ^ 2 Cid 2 C2- 



5F n - nSF 1 

-2, f f lhcomj 

Symmetrising between (1 and (2 we then find 

2 

2-x I fit/- \ \1-x 



if? - usf? = tfdur- / / o/'toi -* ran / (P6<f6 

l/(Ci)-/(c 2 )i>(^) ICi - C2I 




(9) 



6 



The integral is now finite as e — » and the cutoff can be removed, at least for x close to 2. 
Also, notice that everywhere /'(Ci) an d are regular, the integrand vanishes as x — > 2. 

However this is not the case close to the conical singularities. For example, near ( = 0, 
\f'(0\ 2 ~ x ~ |^|( n - 1 )( 2 - x ) ) a nd the limits £ — > and a; — >• 2 do not commute. However 
we would obtain the same limit as x — > 2 in the integral if we restricted the integration 
region to, say, (|G| < p\, IC2I < p) U (|Ci| > — 1 1 , I C2 1 > P _1 ) for any < p < 1. In 
particular, we can take p arbitrarily small, in which case we can accurately approximate 
/'(C) by its asymptotic forms n( n ~ 1 and n(~ n ~ 1 , up to terms which vanish as x — > 2. 
The consequences is that, as x — > 2, we can replace the integral in (|H|) by the analytic 
continuation to x ~ 2 of 

r r \ r ,/-ri-l\2-x\ r ,/-n—l\2—x 

d 2 ( lC l 2 ( 2 . 

This, apart from the £ 4 ~ 2x prefactor, is precisely twice what we would obtain for a single 
branch point. This integral can be evaluated explicitly. Rescaling ( 2 — ^Ci we have 

2n 4 - 2x I(n,x) f \Ci\~ 2 ~ 2n[x ~ 2) d 2 Ci, (10) 

where^l 

T(l + (n + l)(x - 2)/2)r(l - (n - l)(x - 2)/2)r(l - x) 






2— a:| 




2-x 


ICi- 


C2I 


2x 



7T- 



(n + l)(x - 2)/2)r((n - l)(x - 2)/2)T(x) 



Although the equivalence of this result to ([9]) is valid only as x — > 2, we note that it does 
exhibit the required poles at x = 1 and x = 2n/(n — 1), correspond to the short- distance 
divergences already discussed above. However, as x — > 2 we have 

I(n,x) ~ -(7r/4)(n 2 - l)(x - 2) + 0((x - 2) 2 ) . 

The integral over G in f lTU]) . gives, after imposing a short- distance cutoff |£| > e 1 /™, a 
factor 



2vr(e 1 /") 



l/n\ 2n ( x - 2 ) 



71 



+ 0(1) 



2n(x - 2) n(x - 2) 
Putting all these factors together we then find that, for x ~ 2 

r 2 

(i "■ \ 4 



F n -nF 1 = -- (n - n- 1 ) \og(£ / e) + g 2 (n - n- 1 ) (£/ e) 4 ~ 2x l — + 0(x -2) ) +0(g 3 ). (11) 



2 This also shows that the argument also holds in presence of a boundary when the form |Ci — (2\~ 2x 
of the two-point function holds only at short distances. 

3 This follows using the representation | w| 2a = T(— a) -1 J °° u~ 1 ~ a e~ uw w du for each factor and first 
performing the gaussian integral over w. The result can then be reduced to a beta-function integral. 
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For x = 2 this gives an (apparently) uninteresting constant contribution to the Renyi 
entropy, and it is therefore necessary to go the next order, involving an integral over the 
3-point function of the form 

J IC2 - Csl lCs - C1I IC1 - C2I 

where b is the universal coefficient in the OPE 

$(o-$(o) = icr 2K (i+&icr$(o) + •••)• 

Once again, it can be shown that, after subtracting the short-distance singularities as 
Cj ~~ ^ (k, the measure is concentrated on the conical singularities as x — > 2, so the 
result for the interval is twice that found by replacing f(() by ( n . The integral can then 
be performed explicitly. (A similar computation was carried in Refs. [28l [29] for the 
corrections to the free energy of a cylinder, which corresponds to the limit n — > of the 
present calculation.) However, it turns out that the coefficient of the 0(g 3 ) term in (ITTj) is 
determined from the 0(g 2 ) term by Zamolodchikov's c-theorem[23]. To use this, however, 
it is more correct to consider the logarithmic derivative of the free energy with respect to 
the cut-off, which, from ([3]), is proportional to the integral of the trace (O): we see from 
ffTTj) that this takes the form — (c e g(<?)/6)(n — n^ 1 ) where 

CcffG?) = c - 3tt 2 (2 - x)g 2 + 0(g 3 ) . (12) 

The c-theorem[23] states that there exists a function C(g) which decreases along RG flows 
and is stationary at fixed points where it equals the conformal anomaly number of the 
corresponding CFT. Zamolodchikov's analysis also implies that C'(g) oc (l + 0(g 2 ))(3(g), 
where (3(g) = —e(dg/de)\ R , keeping the renormalised coupling fixed. For a perturbed 
CFT, the first two terms in (3(g) are universal[23j [30] : 

-P( g ) = (2-x)g-nbg 2 + 0(g 3 ). 

Thus, up to and including terms 0(g 3 ), all candidates for an interpolating function C(g) 
must agree, and in particular cl ff (g) oc (3(g) to this order. This fixes the coefficient of the 
0(g 3 ) term in flTJ) to be 2tt 3 &H 

However, the result in (JT2"]) disguises the fact that in the neglected higher order terms 
logarithmic dependences on t appear, corresponding to poles at x = 2. These can, 
however, be absorbed by "RG-improving" the expansion, that is, replacing g by g(£) 
where 

tdg(t)ldt=-(3(g(t)). 

If x > 2, that is, the perturbation is slightly irrelevant, then asymptotically g(£) oc f - ^ -2 ) 
and c e ff(£) — > c, with a power- law correction, consistent with our earlier analysis. If x < 2, 

This result is in fact universal in all physical manifestations of perturbative corrections to c, for 
example the conformal anomaly (0) = —cR/\2 in a curved background of scalar curvature R (31 j - 



8 



that is the perturbation is relevant, then, depending on the sign of g/b, either g(£) — > g* 
and c e s(£) — > c new = c — (2 — x) 3 /b 2 + 0((2 — x) 4 ),[29j or g(£) grows beyond the range of 
this perturbative treatment. 

The interesting case is when x = 2, that is, the perturbation is marginal. Then if 
g/b < 0, the perturbation is marginally relevant and g(£) again grows, but if g/b > it is 
marginally irrelevant and g(£) flows to zero, albeit logarithmically slowly. In fact 



* w 1 + nbg \og(£/e) 7rMog(£/e) ' 
Substituting this into ffT2]) we then see that 

= c + b 2 (hg l /e) y + 0(M£/e))- 4 ) ■ (13) 

If we now integrate with respect to e to find the free energy F n — nFi, we find the result in 
02]). Note that this asymptotic value is reached from below, in apparent contradiction to 
Zamolodchikov's c-theorem. However, the correct definition of C(g) in this case is through 
the logarithmic derivative of the entanglement entropy [32], which is ultraviolet finite at 
the fixed points, and in this quantity c is approached from above. We also remark that, 
for fitting any finite £ data, it is better to replace the asymptotic (log(£/e)) 3 in ( {TBI by 



5. Comparison with other studies and discussions 

We have shown that the block entanglement in an infinite Id system (measured by the 
Renyi entropies S^) generically displays standard corrections to scaling of the form 
£-2(x-2) w j^j 1 x > 2 and unusual ones of the form £~ 2x / n (and also the 'combination' 
£2-x-x/n^ w ] 1 i c ] 1 however is never the leading one). We call these corrections unusual 
because of the explicit n dependence of the exponent, a property that seems to contrast 
with RG finite size scaling theory [21] . Clearly there is no contradiction: we showed that 
these terms arise from the conical singularities needed to describe Sjf in a path integral 
formulation. The existence of such geometry-dependent exponents was in fact first noticed 
a long time ago in Ref. [33]. The most surprising effect here is that the unusual correction 
£-2x/n can k e p resen £ a } so f or relevant operators with < x < 2, occasioned by a local 
breaking of scale invariance at the conical singularity. Such effects have probably not be 
seen in earlier studies of corner critical behaviour [34J because they focussed on effective 
values of n < 1 where the exponent of these corrections is larger. 

These unusual effects clearly need direct confirmation from lattice computations. 
Large corrections to scaling have been observed in several numerical studies quoted before, 
but a quantitative analysis has become available only recently [18]. It has been shown 
analytically for the Ising and XX universality class that corrections to the scaling are of 
the form £~ 2 / n [IS], which agrees with our general formula when x — 1. For the Ising 
model, x = 1 corresponds to the energy density operator that indeed we argued to be 
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generated by the conical singularity. For anisotropic Heisenberg chains, the corrections 
to the scaling have been found numerically to be of the form £~~ 2K / n , where K is the Lut- 
tinger liquid exponent, i.e. the most relevant present in the continuum theory. Again this 
perfectly agrees with our result. Similar unusual (i.e. n dependent) corrections have been 
also found in other entanglement measures [351 EE] , but their quantitative understanding 
in the framework of quantum field theory requires further investigation. 

In the case of systems with boundaries, we showed that these unusual corrections are 
of the form i~ x / n . Evidence of such power laws has been reported for n — 1 [191 120] 
both for XX and Heisenberg chains, but a quantitative study for general n is still lacking. 
Preliminary results show that our results are correct [TBI [37] . 

However, our theory of the origin of the £~ 2x / n corrections with x < 2 (and £~ x / n in 
the semi-infinite case) has a definite prediction, so far untested in analytic and numerical 
studies of particular systems: if the origin of these terms is indeed in the local deviation 
from criticality near the conical singularities, then a modification of the lattice action close 
to the singularities should have the effect of changing the amplitudes of these corrections. 
In fact, by tuning the local couplings to a particular value (which may however depend 
on n) it should be possible to eliminate these correction terms altogether. It should be 
stressed that this modification needs to be done locally in space-time, not simply by 
modifying the interaction strengths in the hamiltonian near the ends of the interval. This 
could be carried out by starting with the ground state |0) of the full hamiltonian H, and 
evolving this with the operator exp(— tH'), where H' is H with modified interactions 
in a region 0(e <C £) around the ends of the interval. Here r ~ e/v, with v being the 
coefficient in the quasiparticle dispersion relation u ~ v\k\. The prediction is that the 
leading behaviour of the Renyi entropies measured in the modified state is that same as 
that of the ground state, but that the amplitudes of the 0(£~ 2x ^ n ) correction terms should 
be different. 

It is also worth commenting on the other correction £ 2 - x ~ x / n ^ that is never the most 
relevant one and could be obscured by the others in numerical analysis. For n ~ n c its 
effect should be more important, but an accurate quantitative analysis is difficult because 
the various corrections get too close to each other. However in Ref. [18J it has been 
noticed that for n close to 1 and for A < —0.5 (the anisotropy parameter of the Heisenberg 
chain), the single correction £~ 2K / n does not describe numerical data accurately. This is a 
confirmation of the presence of other corrections to the scaling that could be of the form 
g2-x-x/n ^ u t could also be of very different origin). 

We stress that these corrections in some lattice systems (those described by Luttinger 
liquid theory) show a strong oscillating behaviour [T51 [T^l [T31 EBJ. While the physical 
origin of this effect is due to strong tendency to antiferromagnetic order [19], a general 
proof of this universal form from the continuum Luttinger liquid field theory is still lacking. 

Marginal perturbations deserve separate discussion, because of the logarithmic correc- 
tions to the scaling. It is a well-known (and obvious) fact that log-corrections are hard 
to detect in numerical studies, and in fact, up to now no direct evidence for them has 
been still provided. However, for the isotropic Heisenberg antiferromagnet (that has a 
marginal operator) the relations found in [19] between entanglement entropy (n = 1) and 
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energy of the ground-state (known to have similar kind of log-corrections) is an indirect 
evidence of the correctness of our prediction (cf. Eq. ({TBI ). 

Finally we want to comment on the case of more than one interval, corresponding to 
more than two branch points in the n-sheeted Riemann surface. In this case, already 
the leading term (the analogue of Eq. (TjQ)) is rather involved [391 EES HH H2], but it 
is calculable for integer n for the simplest CFT's [iQl HJJ, [43]. Even in this case, the 
large corrections to the scaling prevent a direct simple analysis [401 HI] and the check of 
the complicated asymptotic forms. The exact knowledge of the correction to the scaling 
exponents, would greatly simplify the analysis. Up to now, it has been shown numerically 
for the Ising model, with periodic boundary conditions, and for n = 2 that the exponent 
is 1/2 [13], compatible with a form such as 2x/n (x = 1/2 and n = 2). Exploiting the 
exact solvability of the model, an accurate analysis [H] shows also for other low values 
of n = 3,4 corrections compatible with 1/n. In the Ising model, the non-local fermion 
operator has dimension 1/2 and could be the origin of such scaling. In fact, to build the 
reduced density of the spin degrees of freedom in the Ising model, one should introduce 
the non local Jordan- Wigner string between the two intervals [431 H21 [44], that plays 
the rule of the Ising fermion. This does not enter in the single interval reduced density 
matrix that is the same for spin and fermion degrees of freedom. As a confirmation of this 
behaviour, we mention that for fermions in the Ising model we have the same correction 
as for the single interval 2/n [44], because the Jordan- Wigner string is not present. For 
the XX model, the same analysis [44] shows unambiguously that the exponent is exactly 
2/n. Thus all evidences confirm that also for two intervals the correction to the scaling 
exponent is of the form 2x/n, even if the example of the Ising model shows that some 
care is needed to fix the appropriate value of x. 
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